GAUSSIAN UPPER DENSITY ESTIMATES FOR SPATIALLY 

HOMOGENEOUS SPDES 
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Abstract. We consider a general class of SPDEs in M. d driven by a Gaussian spatially 
homogeneous noise which is white in time. We provide sufficient conditions on the 
coefficients and the spectral measure associated to the noise ensuring that the density 
of the corresponding mild solution admits an upper estimate of Gaussian type. The 
proof is based on the formula for the density arising from the integration-by-parts 
formula of the Malliavin calculus. Our result applies to the stochastic heat equation 
with any space dimension and the stochastic wave equation with d £ {f , 2, 3}. In these 
particular cases, the condition on the spectral measure turns out to be optimal. 



I. Introduction 

We are interested in establishing Gaussian type upper estimates for the density of the 
mild solution of the following class of SPDEs: 

£u(t,x) = b(u(t,x)) + a{u(t,x))W(t,x), (t,x) G [0,T] x R d , (1.1) 

where T > is some fixed time horizon and £ denotes a general second order partial 
differential operator with constant coefficients, with appropriate initial conditions. The 
coefficients a and b are real-valued functions and W(t,x) is the formal notation for a 
Gaussian random perturbation which is white in time and has some spatially homo- 
geneous correlation (see Section [2.11 for a precise definition of this noise). The typical 
examples of operator £ to which our result applies are the heat operator for any spatial 
dimension d > 1 and the wave operator with d G {1, 2, 3}. 

If £ is first order in time, such as the heat operator £ = — A, where A denotes the 
Laplacian operator on M. d , then we impose initial conditions of the form 

u{0,x) = u {x) xeR d , (1.2) 

for some Borel function u : M. d -> R. If £ is second order in time, such as the wave 
operator £ = JU — A, then we have to impose two initial conditions: 

u(0,x) = u (x), —(0,x)=v (x), xeR d , (1.3) 

for some Borel functions Uq, vq : M. d — > K.. 

The above class of SPDEs has been widely studied in the last two decades. Precisely, 
results on existence and uniqueness of solution in such a general setting have been 
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established in [2H1 EJ QUI E] , while the particular cases of heat and wave equations have 
been studied using several frameworks in [351 EJ |8j [201 HHJ USB EH H] . 

A fruitful line of research developed in some of the above-cited references has been to 
apply techniques of Mallliavin calculus in order to deduce some interesting properties of 
the probability law of the solution at any (t, x) G (0, T] x M. d . In fact, there is a whole 
bunch of results on existence and smoothness of the density for the stochastic heat and 
wave equations, for which we refer to [21 [23 [Tj [201 EH ED E21 EH HE]- Moreover, in 
the paper [23], existence and smoothness of density for the class of SPDEs (II. ip have 
been analyzed, unifying and improving some of the results cited so far. It is also worth 
mentioning that other kind of methods beyond Malliavin calculus can be used to prove 
the absolute continuity of the law of the solution in some particular SPDEs (see e.g. 

Once the existence (and possibly smoothness) of the density of the solution to (II. ip 
is established, one usually gathers at some nice estimates for this density, such as lower 
and upper Gaussian type bounds. Exploiting again techniques of Malliavin calculus, 
this problem has been recently addressed by several authors. Precisely, as far as SPDEs 
with additive noise is concerned, using an explicit formula for the density proved in [21] , 
the main result in [21] says the following (see Thm. 7 therein, and also [25] for related 
results). Let \i be spectral measure associated with the spatial correlation of W, and 
consider the following assumption, which is necessary and sufficient for the existence 
and uniqueness of mild solution to (11. ip (see e.g. [7]). 

Hypothesis 1.1. Let T be the fundamental solution associated to the operator C For 
all t > 0, F(t) defines a non-negative distribution with rapid decrease such that 

$(T):= / / |.7T(£)(£)| 2 n{d£)db < +oo. 

JO JR d 

Moreover, T is a non-negative measure of the form T(t, dx)dt such that, for all T > 0, 

sup T(t,R d ) < +oo. 

0<t<T 



Then, under the above hypothesis, with vanishing initial data, a = 1 and assuming 
that b G C 1 has a bounded derivative, [2H Thm. 7] states that, for small enough t and 
any x G K d , the density p ttX of u(t, x) satisfies, for almost every zGl, 

E \u(t, x) — ml f (z — m) 2 1 E\u(t,x) — m\ f (z — m) 2 

cm exp r^¥(ir I - V,M - exp r^r 

(1.4) 

where m = E(u(t,x)), for some positive constants C\,C2- Note that the term $(£) is 
precisely the variance of the stochastic convolution in the mild form of (11. ip when a = 1. 
As a consequence, this result applies to the stochastic heat equation for any d > 1 and 
the stochastic wave equation in the case d G {1, 2, 3} provided that (see [3 Sec. 3]) 

' /i(d£)<+oo. (1.5) 



On the other hand, in the multiplicative noise setting, such kind of density estimates, 
particularly the lower one, become more difficult to obtain and Nourdin-Viens' density 
formula cannot be applied. This has been already illustrated by Kohatsu-Higa in [16] 
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where, by means of Malliavin calculus techniques, a new method to obtain Gaussian 
lower bounds for general functionals of the Wiener sheet has been obtained. In the same 
paper, the author has applied this result to a stochastic heat equation in [0, 1] and driven 
by the space-time white noise. 

In order to deal with SPDEs beyond the one- dimensional setting, in [26] Kohatsu- 
Higa's general method has been extended to the Gaussian space associated to the un- 
derlying spatially homogeneous noise W. This allowed us to end up with the following 
density estimates for the stochastic heat equation in any space dimension d > 1. Assume 
that b,a G C°° are bounded together with all their derivatives, \&(z)\ > c > for all 
z G R, and for some r\ G (0, 1), it holds 



Then, for all (t, x) G (0, T] x IR d , the density p t)X of u(t, x) verifies, for all z6l, 



where I (t, x) = (T(t) * u )(x), u being the initial data. In the case b = 0, the constant 
C3 would vanish. Note that here, in comparison to (11.41) . the estimates are valid for any 



In fact, let us point out that the upper bound in (11. 7p is much easier to obtain than 
the lower one, and the former comes from the expression for the density popping up 
from the integration-by-parts formula in the Malliavin calculus framework. 

Extending the lower estimate in (11.71) to the general class of SPDEs (II. ip seems to 
be an open problem, for the success in the application of the general strategy of [2E] 
is closely tied to the parabolic structure of the heat equation. However, a much more 
humble objective, which is the one we plan to gather in the present paper, is to tackle 
the upper bound. In fact, we are going to seek the minimal conditions on either the 
coefficients b and a and the spectral measure \i implying that the upper estimate in (11.71) 
remains valid for the general class of SPDEs (11.11) . In particular, we will only need b and 
a to be of class C 2 (and bounded with bounded derivatives) and, for the particular case 
of the heat (resp. wave) equation with any d > 1 (resp. d G {1,2,3}), the condition on 
jj, will be simply (jl.5p rather that (11.61) . More precisely, the main result of the paper is 
the following. We use the notation I (t, x) to denote the contribution of the initial data 
(see (13.21) and (13.31) for its explicit expression in the case of heat and wave equations) 
and suppose that the forthcoming Hypothesis 13.21 is satisfied. 

Theorem 1.2. Assume that Hypothesis \1.1\ is satisfied, and that b,a G C 2 are bounded, 
have bounded derivatives, and \o~(z)\ > c > for all z G R. Moreover, suppose that, for 
some 7 > 0, it holds 




(1.6) 




(1.7) 



T > 0. 





(1.8) 
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Then, for any (t,x) G (0,T] x W 1 , the solution u(t,x) of U.l\) has a density p ttX which 
is a continuous function and satisfies, for all z6R, 



where the constant C3 vanishes whenever b does. 

We remark that, though the above bound does not look exactly Gaussian, it does 
in an asymptotic point of view, namely whenever T is small or z is large. On the 
other hand, we note that, under (11.51) . condition (II. 8p is satisfied for the heat and wave 
equations with 7 = 1 and 7 = 3, respectively (see e.g. [191 Lem. 3.1] and [3H App. A]). 
Similarly, one can also check that the above theorem applies to the stochastic damped 
wave equation with any space dimension (see Example 7 in [TJ Sec. 3]), where condition 
(Oj) is fulfilled with 7 = 3. 

We also point out that our result is not applicable to the case o{z) = z (this would 
be related, e.g., to the parabolic Anderson problem [3]). In fact, in such a case there are 
even very few results on absolute continuity of the law of solutions to SPDEs (see [27]). 
Nevertheless, in the recent paper [13], the authors prove existence and smoothness of the 
density for a stochastic heat equation with a nonlinear multiplicative noise which is white 
in time and with some spatial correlation (much more regular than the one considered 
in the present paper), and with a non-degeneracy condition on the diffusion coefficient 
of the form a(uo(xo)) 7^ for some x G M. d . Their proof is based on a Feynman-Kac 
formula for the solution of the underlying equation. This technique has also been applied 
in [13] to study the density for a stochastic heat equation with a linear multiplicative 
fractional Brownian sheet. 

As mentioned before, the proof of Theorem 11.21 will be based on the expression for the 
density arising from the application of the integration-by-parts formula (see [221 Prop. 
2.1.1]). We point out that this is a well-known method that has been used in other 
contexts (see e.g. [121 E])- As far as the technical obstacles is concerned, the main two 
ingredients needed in the proof of Theorem 11.21 are the following: 

(i) A suitable estimate, in terms of $(£), of the norm of the iterated Malliavin 
derivative in a small time interval (see Lemma 14.11 for details). This will be a 
consequence of a kind of analogous result for the case of the stochastic heat equa- 
tion (see [2S1 Lem. 3.4]) and a mollifying procedure thanks to an approximation 
of the identity which will let us smooth the fundamental solution T(t). 

(ii) A precise control of the negative moments of the norm of the Malliavin derivative 
of the solution, again in terms of $(t) (see Proposition I4.2p . For this, we will 
adapt the proof of [231 Thm. 6.2] to our setting, where the latter allowed the 
authors of that paper to establish that the underlying density is a smooth function 
(under much more regularity on the coefficients though). 

The content of the paper is organized as follows. In Section [21 we rigorously describe 
the Gaussian spatially homogeneous noise W considered in equation (ll.ip . we introduce 
the corresponding Gaussian setting associated to it, together with the main notations 
of the Malliavin calculus machinery. Section [3] will be devoted to recall the definition 
of mild solution to our SPDE (11.11) and summarize the main results on existence and 
uniqueness of solution, Malliavin differentiability and existence and smoothness of the 




(1.9) 
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density. The steps (i) and (ii) detailed above will be tackled in Section HJ Finally, we 
will prove Theorem 11.21 in Section [5j 

With a slight (but harmless) abuse of notation, as already done in this Introduction, 
the notation | • | shall denote either the modulus and norm in M. d . Unless otherwise 
stated, any constant c or C appearing in our computations below is understood as a 
generic constant which might change from line to line without further mention. 



2. Preliminaries 

2.1. Spatially homogeneous noise. Let us explicitly describe here our spatially ho- 
mogeneous noise (see e.g. [7]). Precisely, on a complete probability space (f^J 7 , P), 
this is given by a family W = {W ((/?), <p G C^°(M+ x Mr)} of zero mean Gaussian ran- 
dom variables, where C^°(IR + x M. d ) denotes the space of smooth functions with compact 
support, with the following covariance structure: 

E(W(¥>)W(V0) = / / (v(t,*)*#(t,*))(x)A(dx)dt. (2.1) 
Jo JR d v J 

In this expression, A denotes a non-negative and non-negative definite tempered measure 
on M d , * stands for the convolution product, the symbol * denotes the spatial variable 
and if)(t, x) := ip(t, —x). 

In the above setting, a well-known result of harmonic analysis (see [Ml Chap. VII, 
Theoreme XVII]) implies that A has to be the Fourier transform of a non- negative 
tempered measure \i on M. d , where the latter is usually called the spectral measure of the 
noise W. We recall that, in particular, for some integer m > 1 it holds 

/ t, — L.o\ M^O < +°° 

j Rd (i + iei 2 ) mPV 

and, by definition of the Fourier transform in the space iS'(IR d ) of tempered distributions, 
A = T[L means that, for all <fi belonging to the space iS(M d ) of rapidly decreasing C°° 
functions, 

/ <f>(x)A(dx) = [ 

Therefore, we have 

eN 2 ) = / / \^m)\ 2 m)dt. 

JO JR d 

A typical example of space correlation is given by A(dx) = f(x)dx, where / is a non- 
negative function which is assumed to be integrable around the origin. In this case, the 
covariance functional (12.11) reads 




(p(t, x)f(x — y)ip(t, y) dydxdt. 



The space-time white noise would correspond to the case where / is the Dirac delta at 
the origin. 
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2.2. Gaussian setting and Malliavin calculus. We are going to describe the Gauss- 
ian framework which can be naturally associated to our noise W and introduce the 
notations involved in the Malliavin calculus techniques. 

To start with, let us denote by H be the completion of the Schwartz space <S(IR d ) 
endowed with the semi-inner product 

((/> 1 *<p 2 )(x)A(dx)= [ FMOFMSKdZ), 1 ,0 2 G«S(R d ). 
Jm. d Jm d 

As proved in [7J Example 6], we remind that the Hilbert space TL may contain distribu- 
tions. 

Let T > be a fixed real number and define %t '■= L 2 ([0, T]; %). Using an approx- 
imation argument, our noise W can be extended to a family of mean zero Gaussian 
random variables indexed by %t (see e.g. [6] Lemma 2.4]). With an innocuous abuse of 
notation, this family will be still denoted by W = {W(g),g G Ht}- Moreover, it holds 
^{W(gi)W(g2}) = (gi,g2)u T i f° r au 9\i9i e Ht- Thus, this family defines an isonormal 
Gaussian process on the Hilbert space Tir and we shall use the differential Malliavin 
calculus based on it (see e.g. [221 133])- 

As usual, we denote the Malliavin derivative operator by D. Recall that it is a closed 
and unbounded operator defined in L 2 (Q) and taking values in L 2 (Q; Ht), whose domain 
is denoted by D 1 ' 2 . More general, for any integer m > 1 and any p > 2, the domain of 
the iterated Malliavin derivative D m in L P (Q) will be denoted by D m ' p , where we remind 
that D m takes values in W^Uf* 1 ). We also set ©°° = H p >i H meN B m 'P. The space 
D m,p can also be seen as the completion of the set of smooth junctionals with respect to 
the semi-norm 

{m ^ p 

E(iFn+E E (ii^ F &r)) • 

For any differentiable random variable F and any r = (r l5 ...,r m ) G [0,T] m , D m F(r) is 
an element of %® m which will be denoted by D™F. 

A random variable F is said to be smooth if it belongs to D°°, and a smooth random 
variable F is said to be non- degenerate if G n p >iL p (tt). Owing to [22| Theorem 

2.1.4], we know that a non-degenerate random variable has a C°° density. 

For any t G [0, T], let T t be the a- field generated by the random variables {W s (h), h G 
U, < s < t} and the P-null sets, where W t {h) := W{l m h). 



3. Spatially homogeneous SPDEs 

We gather here a general result on existence and uniqueness of mild solution for our 
SPDE f 1 1.11) and the main results on Malliavin calculus applied to it, namely Malliavin 
differentiability and existence and smoothness of density. As usual, we will also focus on 
the main examples of application that we have in mind, which are the stochastic heat 
and wave equations with d > 1 and d G {1, 2, 3}, respectively. 

We recall that, by definition, a mild solution of (11. ip is an T t - adapted random field 
{u(t, x), (t,x) G [0,T] x R d } such that the following stochastic integral equation is 
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satisfied: 

u(t, x) — Io(t, x) + / / T(t — s, x — y)a(u(s, y)) W(ds, dy) 

JO JR d 

+ [ [ b(u(t - s,x - y))T(s,dy)ds, P-a.s., (3.1) 

Jo Jl$L d 

for all (t, x) G [0, T] x M. d . Here, T denotes the fundamental solution associated to C and 
Io(t, x) is the contribution of the initial conditions, which we define below. 

The (real-valued) stochastic integral on the right-hand side of (13. ip is understood with 
respect to the cylindrical Wiener process that can be naturally associated to our spatially 
homogeneous noise W (see [23J E] and also |3SJ E]). In particular, we will assume that 
Hypothesis II. II is satisfied. Concerning the last integral on the right-hand side of ( 13. ip . 
we point out that we use the notation 'T(s, dy)" because we will assume that T(s) is a 
measure on Mr. 

As far as the term I (t,x) is concerned, if £ is a parabolic-type operator and we 
consider the initial condition (11.21). then 



I {t,x) = (T{t)*u ){x) = u {x-y)T{t,dy). (3.2) 

Jw. d 

On the other hand, in the case where C is second order in time with initial values (11.31). 



I (t, x) = (T(t) * v ) (x) + j t (T(t) * u ) (x). (3.3) 

Example 3.1. Owing to the considerations in [7J Section 3] (see also [231 Examples 4.2 
and 4.3]), in the case of the stochastic heat equation in any space dimension d > 1 and 
the stochastic wave equation in dimensions d = 1,2,3, the fundamental solutions are 
well-known and the conditions in Hypothesis 11.11 are satisfied if and only if 

1 ■ n(dg) < +oo. (3.4) 



1 + 



We shall consider the following assumption on the initial conditions. In the case of the 
stochastic heat equation in any space dimension (resp. wave equation with dimension 
d — 1, 2, 3), sufficient conditions on u (resp. uq, Vq) implying that the hypothesis below 
is fulfilled are provided in [61 Lemma 4.2]. 

Hypothesis 3.2. (t,x) !->■ I (t,x) is continuous and sup( t)a .) e [ 0)T ] xR d |/ Q (t, x)| < +oo. 

The following well-posedness result, which is a quotation of [61 Thm. 4.3], is a slight 
extension of the results in [7j. 

Theorem 3.3. Assume that Hypotheses \1.1\ and \3.2\ are satisfied and that a and b are 
Lipschitz functions. Then there exists a unique solution {u(t,x), (t,x) G [0,T] x IR d } of 
equation \3. Moreover, for all p > I, 



sup E{\u{t,x)\ p ) < +oo. 

{t,x)e[0,T]xR d 
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Let us now deal with the Malliavin differentiability of the solution u(t, x) of (13. ip . For 
this, we consider the Gaussian context described in Section 12.21 The following propo- 
sition summarizes a series of results in [T9(, [321 E3]. For the statement, we will use the 
following notation: for any m G N, set s := (si, . . . , s m ) G [0,T] m , 
(R d ) m , s(i) : = (si, . . . , Sj_i, Si+i, . . . , s m ) (resp. z(i)), and, for any function / and vari- 
able X for which it makes sense, set 

A m (f,X) := D m f(X) - f'(X)D m X. 

Note that A m (f,X) = for m — 1 and, if m > 1, it only involves iterated Malliavin 
derivatives up to order m — 1. 

Proposition 3.4. Assume that Hypothesis \1.1\ is satisfied and, for some m G NU {oo} ? 
cr, b G C m (R) and their derivatives of order greater than or equal to one are bounded, 
Then, for all (t, x) G [0, T] x ~R d , the random variable u(t, x) belongs to W' p for any 
j = 1, . . . , m and p > 1. Furthermore, for any j G {1, . . . , m} and p > 1, the iterated 
Malliavin derivative D^u{t,x) satisfies the following equation in L P {VL\ T-i^): 

D 3 u(t,x) = Z j (t,x) 

+ / / T(t ~ s,x - y^i^uis^)) + D^u(s,y)a'(u(s,y))]W(ds,dy) 

JO JR d 

+ j / [A J (6, u(t — s, x — y)) + D^u{t — s, x — y)b'(u(t — s, x — y))] T(s, dy)ds, 

(3.5) 

where Z^(t,x) is the element of L p (fl;'H'^ 3 ) given by 

j 

Z J {t,x)s,s = z^2 T ( t - Si,x- dZi)Dl^-^a(u(si, Zi)). 
i=i 

A detailed description of the construction of Hilbert-space-valued stochastic integrals 
as the one in equation (13.51) can be found in [231 Sec. 3]. Indeed, as proved in [HI Sec. 
3.6], these kind of integrals turn out to be equivalent to Hilbert-space-valued stochastic 
integrals a la Da Prato and Zabczyk [5]. 

The above Proposition 13.41 can be used to obtain the following results on existence 
and smoothness of the density for the solution u(t,x). They are direct consequences of 
Theorems 5.2 and 6.2 in [23], with the only difference that the latter consider vanishing 
initial conditions. 

Theorem 3.5. Assume that Hypotheses \1.1\ and \3.S\ are satisfied, b,o~ G C 1 have a 
bounded derivative, and \&(z)\ > c > for all z G R. Then, for all (t,x) G (0, T] x R d , 
the random variable u(t, x) has a law which is absolutely continuous with respect to the 
Lebesgue measure. 

Theorem 3.6. Assume that Hypotheses \1.1\ and \3.S\ are satisfied, a, b G C°° and their 
derivatives of order greater than or equal to one are bounded, and that \o~(z)\ > c > 
for all z G R. Moreover, suppose that, for some 7 > 0, 

CF<[[ |JT(s)(0| 2 ri<%)ds, * 6(0,1). (3.6) 

JO JR d 
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Then, for every (t,x) G (0, T] x IR d , the law of the random variable u(t,x) has a C°° 
density. 

As commented in the Introduction, both results apply to the stochastic heat equation 
with d > 1 and the stochastic wave equation with d G {1,2,3} provided that (13.41) 
is satisfied, since condition (13. 6p holds for these examples with 7 = 1 and 7 = 3, 
respectively. 

4. Auxiliary results 

This section is devoted to prove the main two ingredients needed in the proof of 
Theorem II .21 The first one establishes a suitable uniform bound for the norm of the 
iterated Malliavin derivative of the solution u(t, x) in small time intervals. The second 
one deals with the negative moments of the corresponding Malliavin matrix, which here 
simply reduces to the norm of the Malliavin derivative of u(t, x). 

Lemma 4.1. Let < a < e < T and p > 1. Assume that Hypotheses li.il and \3.2\ 

are satisfied and that, for some m G N, the coefficients b, a belong to C m and all their 
derivatives of order greater than or equal to one are bounded. Then, there exists a positive 
constant C, which is independent of a and e, such that, for all 5 G (0, e — a]: 



sup 

(r,y)e[e— S,e] xl 



Ef||I>»"«(r,y)||* ) <Cm JP , (4.1) 

! \ H e-S,e J 

for all j G {1, . . . , m}, where we remind that, for all t > 0, 

JO JR d 

Proof. It is similar to that of j26l Lem. 3.4], where a conditioned version of this result 
for the stochastic heat equation has been proved. Precisely, as already pointed out in 
[26| Rmk. 3.5], in our general setting we need to smooth the fundamental solution T as 
follows. Let if) G C^°(IR d ) be such that ip > 0, its support is contained in the unit ball of 
R d and f Rd i/j(x)dx = 1. For n G N, set ip n (x) := n d ifj(nx) and, for all t, T n (t) := tj) n *r(t). 
It is well-known that T n (t) belongs to 5(IR d ). 

Let us now consider {u n (t,x), (t,x) G [0,T] x R d } the unique solution of 



u n (t,x) = I (t,x) + / T n (t - s,x - y)a(u n (s,y))W(ds,dy) 

JO JR d 

+ 11 b(u n (t - s,x - y))T(s,dy)ds. 
Jo Jm d 

Since T n (t) is a smooth function (such as in the case of the heat equation), we can mimic 
the proof of [26], Lem. 3.4], so that we end up with estimate ( 14. ip with u replaced by 
u n . Indeed, we should remark here that the term involving the pathwise integral with 
respect to T(s, dy)ds does not cause any problem since we only need to use that T(t, M. d ) 
is uniformly bounded in t, which is part of Hypothesis 11.11 

On the other hand, a direct consequence of the proofs of [32, Thm. 1] and [231 Prop. 
6.1] is that, for all (£, x) G [e — 5, e] x M d and j G {1, . . . , m}, 

D j u{t,x) = L 2 (n ; Hfl 5e ) - lim D j u n {t,x). 
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Therefore, writing down the corresponding convergence of norms and taking supremum 
over [e — 5, e] x R d , we conclude the proof. □ 

Proposition 4.2. Assume that Hypotheses li.il and Iff.ifl are satisfied, that b,a are C 1 
functions with bounded derivatives and that \o~(z)\ > c > for all z G R. Moreover, 
suppose that, for some 7 > 0, 

CV<fl IJT0X0I 2 KdOds, te(0,l). (4.2) 

JO JR d 

Then, for any p > 0, there exists a constant C > such that, for all (t, x) G (0, T] x M. d , 

E(\\Du{t,x)\\£)<C*{t)-*. 

Proof. The proof's structure is analogous as that of the proofs of [23j Thm. 6.2] and 
[26"t Prop. 4.3], so we will only sketch the main steps. 

First, owing to [251 Lem. 2.3.1], it suffices to check that, for any q > 2, there exists 
e o — £ o(q) > such that, for all e < e , 

P{Q{t)- 1 \\D«{t,x)\& T <e}<Ce*. (4.3) 

Note that the Malliavin derivative Du(t, x) verifies the following equation in (take 
m = 1 in (133): 



Du(t, x) = cr(u(-, *))r(£ — •, x — *) 
+ 



/ / r(t-s,x- y y(u(s,y))Du(s,y)W(ds,dy) 

J0 jR d 

+ / / b'(u(s,x — y))Du(s,x — y)T(t — s,dy)ds. 



'0 



Then, for any small 5 > one proves that 

P{^ty^Du(t,x)f nT <e}<F{^(t)- 1 I(t,x;5)>c^(t)- 1 ^(5)-£} 

< (c$(t)- 1 $(5) - e)~ P <5>(t)- p E(\I{t,x;5)\ p ), (4.4) 
where I(t,x;S) := WR^t, x; S)\\^ tst + \\R 2 (t,x;6)\\^ tst and 

R x {t,x-8) := I [ r(t-s,x-y)a'(u(s,y))Du(s,y)W{ds,dy), 
J Js. d 

R 2 (t,x;5) := / / b'(u(t — s,x — y))Du(t — s,x — y)T(s,dy)ds. 

J JR d 

Using the above Lemma 14.11 and applying standard integral estimates, one checks that 

E(\I(t,x;S)\ p ) < C<5>(5) p (<5>(5) p + y{5) p ), 

where we have set 

*(a) := / T(r,R d )dr. 
Jo 

Thus, going back to (14.41) . we obtain 
F^ity^Du^x)^ < s} < C (c^it)- 1 ®^) - sy P <S>(t)- p $(5) p ($(6) p + y(6) p ) . 
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At this point, taking a small enough e if necessary, we can choose 6 = 5(e) such that 

U(t)-'m=e. (4.5) 

Hence, we have 

P{$(t)" 1 ||lMt,aO||* T <e} < C($(5) p + M>(<5) p ). 

Note, on the one hand, that condition (14.51) implies $(5) < C$(T)e < Ce. On the other 
hand, by Hypothesis ll.il we have *&(5) < CS. Hence, the assumption (14.21) and what we 

have just said above let us infer that ^(5) < Ce~< . Therefore, 

P{$(t)- 1 ||Dn(t,x)||^ T <4<c(^ + ^), 
so taking p = q(j V 1) we conclude that (14. 3 p is satisfied. □ 



5. Proof of the main result 

In this section, we are going to prove Theorem 11.21 On the one hand, we note that 
Proposition 13.41 implies that, for any (t, x) G (0,T] x M. d , the random variable u(t,x) 
belongs to D 2 ' p for all p > 1. Moreover, an immediate consequence of Proposition 14.21 
is that the Malliavin matrix associated to u(t, x) has negative moments of all orders. 
Thus, applying a general criterion of the Malliavin calculus (see e.g. [2"2"l Prop. 2.1.5] or 
[TT1 Thm. 4.1]), we obtain that the law of u(t, x) has a density and it is a continuous 
function. 

On the other hand, as explained in the Introduction, the proof of (II. 9p is a matter of 
following exactly the same arguments as in [261 Sec. 5] and invoking the results of the 
previous section. Let us sketch the main steps to follow. 

To start with, we use the formula for the density arising from the application of the 
integration-by-parts formula in the Malliavin calculus context (see e.g. [221 Prop. 2.1.1]). 
Precisely, denoting the density of u(t, x) by p t>x , we have 

„ / . / Du(t,x) \\ 

pUv) = E [iMt,* >y} 6 [ llDu ^ x)rHT )), v e R, 

where here 8 denotes the divergence operator or Skorohod integral, that is the adjoint 
of the Malliavin derivative operator (see [2"2"l Ch. 1]). 

Next, taking into account the equation satisfied by u(t,x) (i.e. (13.11) ) and applying 
[221 Prop. 2.1.2], we obtain 

PtM <CP{|M t | > \y-I Q (t,x)\-c 3 T}* 

x |e (\\Du(t, x) H^) + (E\\D 2 u(t, x) ||£ ?2 ) 5 (E\\Du(t, x) ||^f ) J j , 

(5.1) 

where a,(3,q are any positive real numbers satisfying ^ + ^ + ^ = i - 111 the above 
expression, M t denotes the martingale part of the solution u(t, x), that is 

M t = I I T{t - s,x - y)a{u{s,y))W{ds,dy), 

JO JR d 
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and the term c^T comes from the fact that, due to Hypothesis 11.11 and the boundedness 
of b, for all (t,x) G (0,T] x R d , 



/ / r(t- s,x -y)b(u(s,y))dyds 

J0 JR d 



<c 3 T, 



P-a.s. 



In order to estimate the terms in (15. ip . we first apply the exponential martingale 
inequality in order to get a suitable exponential bound of the probability in (15.1 j) (using 
that (M)t < C $(£)), and then we conveniently apply Lemma [4.11 and Proposition 14.21 
Thus 

(\y-I (t,x)\-c 3 T)^ 



Pt, 



cxp 



c 2 $(t) 



y e 



where the constants ci, C2, C3 do not depend on (t, x), so we conclude the proof of Theorem 

o □ 
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